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Abstract. We study the blowup criterion of smooth solution to the Oldroyd models. Let 
(u(t, x), F(t, x) be a smooth solution in [0,T), it is shown that the solution (u(t, x), F(t, x) 
does not appear breakdown until t = T provided Vu(t,x) G L 1 ([0, T] ; L°° (K n )), n = 2, 3. 
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1. Introduction 

In this paper, we consider the blowup criterion of smooth solution to the incompressible Oldroy 
model in the two and three dimensional space: 

!d t u - vAu + u ■ Vu + Vp = V • (FF l ), 
d t F + u-\7F = VuF, 
divu = 0, 

for any t > 0, x E W l , n = 2, 3, where u(t,x) is the velocity field, p is the pressure, fj, is the 
viscosity and F the deformation tensor. We denote (V--F)j = d Xj Fij for a matrix F. The Oldroy 
model (jl.ip describes an incompressible non-Newtonian fluid, which bears the clastic property. 
For the details on this model see [7]. 

The local existence and uniqueness of the Oldroy model on entire space R™ or a periodic 
domain was established by Lin etc. in [7J, where the global existence and uniqueness of smooth 
solution with small initial data was also established see also jS] . The wellposedness on a bounded 
smooth domain with Dirichlet conditions was established by Lin and Zhang in [S]. 

We remark some properties of the deformation tensor. Let x be the Euler coordinate and X 
the Lagrangian coordinate. For a given velocity field u(t,x) the flow map x(t,X) is defined by 
the following ordinary differential equation 

{±x(t,X) = u(t,x(t,X)), 
\x(0,X) = X. 

The deformation tensor is F(t,X) = J^(i,X). In the Eulerian coordinate, the corresponding 
deformation tensor is define as F(t, x(t, X)) = F(t, X). Differentiating its both sides with respect 
to t by chain rule one obtain the second equation of which says that d t Fij + Uk ■ d Xk Fij = 

dx k UiFkj for i, j = 1, 2, • • • , n, in the th entries, where we use the Einstein summation 

convention that the repetition index denotes sum over 1 to n. 

If div -F(0, x) — 0, then from the second equation of Oldroy (jl.ll) we have 

(1.2) 5 t (V-F t )+u-V(V-jP rt ) = 0. 

Therefore, V ■ F* = for any t > 0. 
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Denote the ith column of F as F.i, then V • (FF*) = F. t • V-F.j by the fact V • F* = 0. So the 



system (jl.ip can be rewritten in an equivalent form 
(1.3) 



d t u - vAu + u ■ Vu + Vp = F.j • Vf, 
dtF.fc + it ■ VF. fc = F. fc • Vu, k = 1, • ■ 
divu = 0, divF = 0. 



. n, 



In reference [7j, Lin, Liu and Zhang obtained the local existence and uniqueness of smooth 
solution for smooth initial data, and had a blowup criterion. 

Theorem (Lin, Liu and Zhang) For smooth initial data (u ,F ) E H 2 (R n ), there exists a 
positive time T = T(||u ||#2, ||F ||#2) such that the system (jl.ip possesses a unique smooth 
solution on [0,T] with 

(u,F) G L°°([0,T];H 2 (R n ))nL 2 ([0,T];H 3 {R n )). 
Moreover, if T* is the maximal time of existence, then 



T* 



Vu\\ H - 2 ds = 



In reference [3], Hu and Hynd study the blowup criterion for the ideal viscoelastic flow, which 
is the Oldroy system (jl.lj) in the case of /i = 0. They showed an Beale-Kato-Majda p] type 
blowup criterion that the smooth solution to the Oldroy flow do not develop singularity for t < T 
provided that 



|V x u\\ L , 



|V x F. k \\ L ^, m As < +oo. 



From the modeling of Oldroy system we know that the deformation tensor can be determined 
by the velocity u of the flow. Therefore we consider the blowup criterion of smooth solution by 
means of only HVu^. In fact, Zhao, Guo and Huang constructed a set of finite time blowup 
solution in two dimension case: 

t 



u(t, x) 



F(t,x) 



xifo 



£2/0 



f - 
diag 



g-a 



a- /3 



fot 



p(t,x) 
a + [3 



{axl-Px 2 2 )f§ 
(/3-«)(l-f^f/oi) 2: 



fat 



If ^r§/o > 0, a+P ^ and a — P ^ 0, then the above solution will blow up at time T" 
We see that 

||Vv(t)||oodf = +00. 

There are other types of blowup criteria of smooth solutions to the Oldroy models, for example 
[BJE]- To this end, we state our main results. 

Theorem 1.1. Let u G H 2 (R n ) and F £ H 2 (R n ) with V • u = V • F. kfi = for k = 1, • • • , n. 

Assume the pair (u,F) G L°°([0,T];H 2 {R n )) n L 2 ([0, T];H 3 (R n )) is a ' smooth solution to the 
Oldroy system U.3\) . Then the smooth solution do not appear breakdown until T* > T provided 
that 



_ a—B 



(1.4) 



WVu^Woodt < +00. 



Remark 1.1. For the local smooth solution (u,F) G L°°([0, T];H 2 (R n )) n L 2 ([0, T];H 3 (R n )), if 
T* is its maximum existence time, then ||Vu(t)||oo(if = +00. 
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In the second section we will prove the Theorem II .11 for the case n = 2, which can be done by 
energy estimates. The L 2 and H 1 energy estimates are the same for the case n — 2 and n = 3. In 
the H 2 energy estimate, we use the Sobolev interpolation inequality ||V.F||| < C\\ Vi^|| 2 1| AF|| 2 . 

1 3 

In case n = 3, however, the inequality is ||VF||4 < C||VF||| |AF||| which does not match the 
H 2 energy estimate, because it will result in the appearance of the term |AF||2 that the power is 
higher that the left hand side. We obtain the H 2 energy estimate of u by virtue of the momentum 
equation, combining the H 2 estimate of u and F again with the estimate of llVi 7 ^^ we grasp 
the H 2 energy estimate of u and F finally. The section three will devote to the proof of the case 
n = 3. 

In this paper C denote a harmless constant which may be dependent on dimension n, the 
norm of initial data, the viscosity /i, but not dependent on the estimated quantity. We denote 
the LP norm of a function / by ||/|| p or ||/||lp. We denote the derivative with respect to Xi by 
di or d Xi . We also use ft to denote the derivative of / with respect to t. 

2. Proof of the case n = 2 
(1) L 2 -energy estimate and LP estimate of the deformation tensor F 

The L 2 -energy estimate can be easily obtained by the standard L 2 inner product process. 
l^(\\ u \\l + \\ F -k\\l) + H|Vu||2 = (F-i ■ VF.i,«) + (F. k ■ Vu,F. t ) = 0. 

So we have 

(2.1) \\u\\ 2 + \\F\\ 2 + 2^ /* HV^II^ = ||uo||3 + ||F(0)|| 2 . 

Jo 

Multiplying both sides of the second equation of (jl.3|) by p\F.k\ p ~ 2 F.k for 2 < p < oo and 
integrating both sides on W l it follows that 

(2-2) ^ t \\Fk\\ p p <p\\Vu\U\F\\l. 



Summing up the estimate (|2.2|) with respect to k one has 

(2.3) \\F\\ p <\\F \\ p expic(n) [ ||Vw(s)|| 00 ds 

Let p — > oo, we have 



(2.4) 



l^lloo < ||^o||ooexp|c(n)^ ||Vu(a)||ood« 



(2) ij 1 -energy estimate 

We differentiate the equations (|1.3I) with respect to Xi, then multiply the resulting equations 
by diU and diFj for i = 1, 2, integrate with respect to x and sum them up. It follows that 

5^(IMI + MI) + mII^v u ||1< 

\(d iU ■ Vu,d lU )\ + KdiFk ■ VF. k ,d % u)\ + \(d iU ■ VFj^iF^l + {(diFj ■ Vu^iF.^l 

where use has been made of the facts 

(u ■ VdiU, diu) = (u ■ i), I-.,} = (Vdip, d l u) = 0, 
{F.k ■ \7d t F k ,d lU ) + (F.j ■ V8 t u, diF.j) - 0. 

Noting that 

\(d iU - Vu,diu)\ < || V^Hooll V^|tl, 
\{d t F. k -VF. k ,d t u% {(diU-VF^diF^l, \(d i F. r 'S7u,8iF. j )\ < \\ Vu||oo|| VF||1. 
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So 

^(V«lll + l|VF]|) + »\\D 2 u\\l < C\\Vu\U\\Vu\\l + \\VF\\l). 
Gronwall's inequality implies 

(2.5) ||V^[|1 + \\VF\\l + 2Ai jT ||£> 2 i/[||ds < (UVuolll + ||VF(0)[[|) exp | J* C|| Vu(s) lUdsj . 

(3) iJ 2 -energy estimate 

Applying operator A on both sides of (11.31) . we have 
(2.6) 

jd t Au - uA 2 u + Au ■ V« + u ■ VAu + 2d t u ■ Vd t u + VAp = AF. k ■ VF. k + F k VAF k + 2d l F. k ■ Vd l F. k 
\d t AF. k + Au ■ VF. k + u ■ VAF.fe + 2d l u ■ Vd l F. k = AF. k ■ Vu + F. k ■ VAu + 2d t F. k ■ Vd l u. 

Taking the L 2 inner of equation (|2.6[) with Au and AF. k and summing them up, one can obtain 
that 

~(||Au||§+||AF||!) + M ||AV«||2 
t 2 7 \ < \(Au ■ Vu, Au)| + 2\{diU ■ VdiU, Au)\ + \{AF. k ■ VF k ,Au)\ 

+ 2\(diF k ■ V^F.fe, A«)| + |(Au • VF. fc) AF k )\ + 2\{d lU ■ Vd. t F k , AF k )\ 
+ \(AF k -Vu,AF k )\ + 2\(d i F k -Vd i u,AF k )\. 

Here use has been made of the the facts that 

(u ■ VAu, Au) = 0, (u ■ VAF. fc , AF k ) = 0, 
(F. k • VAF. k , Au) + (F. k ■ VAu, AF. k ) = 0. 

Noting that 

(Am • Vu,Au)\, \(d t u ■ Vd l u,Au)\ < C\\ Vu\\oo\\ Au\\j, 
\(AF k -Vu,AF k )\, \(diU-VdiF. k ,AF. k )\ < C[|Vu||oo||AF||2. 
\(AF k ■ VF. k ,Au) + 2{d t F. k ■ VdiF k ,Au)\ 

= | - (diF k ■ VF k ,diAu) - (diF k ■ VAu,d,F k )\ < C\\VAu\\ 2 \\VF\\l 

< J||VAu||l + C||VF|||||AF||l, 
where we have used the Sobolev interpolation inequality 

||VF||1<C||V^|| 2 ||A J P|| 2 . 
Arguing similarly as the above, one has 

\(Au-VF. k ,AF k )\ = KdiAu-VF^diFk)] < ^\\D 3 u\\j + C\\VF\\ 2 2 \\AF\\l 

o 

\{diF. k ■ Vd t u, AF. k )\ = \{didjF. k ■ Vd 3 F k , d iU )\ + \{d t F. k ■ Vd^u, d,F. k )\ 
< C\\Vu\U\AF\\ 2 + |||VAw||l + C\\VF\\ 2 \\AF\\ 2 . 
Inserting the above estimates into estimate (|2.7j) . it can be derived that 

~(l|Att||l + \\AF\\ 2 2 ) + |||AV«||5 < C||Vu|&(||Au||2 + ||AF||1) + C\\ VF\\ 2 2 \\ AF\\ 2 2 . 

Gronwall's inequality implies that 
(2.8) 

||A«||| + \\AFg + M ^ \\AVu( s )\\ 2 ds < (\\Au g + ||Af (0)||1) exp j^cxp^ Ct||V«||ood«|. 

(4) Higher derivative estimates. 
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Next we derive the higher derivative estimate of u and F . For this purpose we need the 
following commutator estimate. 

Proposition 2.1. (Kato and Ponce [3], [5]) Let 1 < p < oo andO < s. Assume that /, g S W s ' p , 
then there exists a abstract constant C such that 

(2.9) \\{J s J}g\\ P < C(||V/|| P1 ||<?|| w .-i, P2 + H/lk-alMU) 

(2-10) l|[A s ,/].g|| P < CdlV/llpJIA-^Hp, + ||AVIUIslU 

with 1 < P2,P3 < oo such that 

1 _ 1 11 1 

P Pi P2 P3 Pa 

where [A s J}g = A s {fg) - fA s g and A s = (-A)2, J = (1 - A) 1 / 2 . 

Applying A s on both sides of (|1.3p and taking the inner product with A s u and A S F, it can be 
derived that 

~(||A'u||i + ||A«F. fc ||l) + /i||A' +1 u||2< 
( 2 - n ) \{A s (u- Vu) - u- VA s u,A s u)\ + \(A s (F. k ■ VF. k ) - F. k ■ VA s F. k: A s u)\ + 

\(A s (u ■ VF. k ) - u • VA s F k ,A s F k )\ + \(A s (F. k ■ Vu) - F k ■ VA s u, A s F. k )\, 

where we have used the facts 

(F k ■ VA s F. k , A s u) + (F k ■ VA s u, A s F k ) = 0, 

(u ■ VA s F. k ,A s F. k ) = (u • VA S «, A s u) = 0. 
The commutator estimate (|2.10j) implies that 

||A s (u ■ Vu) - u ■ VA s u\\ 2 < UVulloollA'ulla, 
||A s (F. fe • VF k ) - F. k ■ VA s F k \\ 2 < IjVFjUjA^Ua < \\VF\\ B .-i \\A S F\\ 2 , 
\\A s (u ■ VF k ) - uVA s F k \\ < llVulloollA'Flla + ||F|U|A S+1 U || 2 , 

\\A s (F k ■ Vu) - F k VA s u\\ < ||V«||oo||A'F||a + I^UIA^a, 

where the Sobolev embedding H a - 1 (R n ) ^ L°°(R") for s > 1 + § is applied. 
Inserting the above estimates into estimate (|2.1ip . it follows 

(212) ^(l|A^ + ||A^||i) + f||A^||< 

CdlVulloc + ||VF|| a + ||A*u||a + ||F||^)(||A*u||i + ||A'F||2), 

where we have used the fact 

||VF|| ff «-x||A*F|| 2 ||A s U ||a < ||VF|| 2 (||A'F||2 + \\A s u\\l) + ||A S F|| 2 ||A^|| 2 . 

So, for s > 3, applying Gronwall's inequality to (|2 . 12|) . by induction for it's estimate, we obtain 
the higher derivative estimate: 

\\A s u\\l + \\A s F\\l + ^ f \\A s+l u\\l&s < 
Jo 

(\\uo\\h + \\F(0)\\U ex P | jT C(l|V«||oc + l|VF||a + ||A S U || 2 + ||F|| 2 0o )d S |. 
Therefore, we complete the proof of the case n = 2. 
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3. Proof of the case n = 3 

In the three dimensional case the L 2 and H 1 energy estimates are the same as the case of 
dimension two. To estimate the H 2 energy estimate we need the following estimates. 

Multiplying the first equation of (|1.3[) by u t and integrating both sides over R 3 with respect 
to x, and noting div u — 0, it follows 

f^llVull! + IKII3 < |(« • V«,ttt)| + \(R k ■ VF k ,u t )\ 

<5lKlll + ^IHl2ol|v«||l + q|vj'||l||F||2 . 

Integrating both sides with respect to t it yields 

(3.1) H|V«||1+ / \\utgds < mHVwqHI + sup \\Vuf 2 f \\u\\ 2 H2 ds+ f \\FWUVFWlds 

JO 0<s<t Jo Jo 

where the Sobolev embedding i/ 2 (M 3 ) L°°(R 3 ) has been used. 

Differentiating the first equation of (|1 .3[) with respect to t, we arrive at 

(3.2) u tt - ^Au t + u t -Vu + u-Vu t + Vp t = F. kt ■ VF. fe + F. k ■ VF. kt . 

Taking L 2 inner product of the equation (|3.2[) with respect to ut, it can be similarly derived that 

^Wll + Al^Ml < \\Vu\U\u t }\l + 2\\F\U\Vu t h\\Fth 

<fl|v Ut ||I + ||v u |U]K||I + 6-11^11^111. 

Applying the Gronwall' inequality, it yields 

(3.3) IKHI + // [ \\Vu t \\lds< (\\u t (0)\\ 2 + C /Vl| 2 ooll^ll2d S )exp( f WVu^ds 

Jo Jo I Jo 

It need still to estimate H-FtH!- From the second equation of (|1.3[) it can be derived that 
ll^lll < llFtllaNUIlVFHa + ||F || 2 ||F|U|Vu|| 2 



So we arrive at 



< ^11^111 + CH| 2 ocl|VF||2 + C||F||^||V W |||. 



F t \\ 2 , < C\\u\\ 2 \\VF\\i + C\\F\\* llVdl 2 



H\\2 ^ ^ II "Mooll vr || 2 -r^\\r Noo || vu|| 2 
Inserting it to the estimate (|3.3I) we obtain the estimate of ||w*||2: 



(3.4) \\u t \\ 2 + 



H [ ||Vitt||lds < C(t) < oo, 
Jo 



where C(t) is explicit increasing function of i dependent on J Q ||Vu||oods. From the first equation 
of (|1.3p . Vp can be solved by Riesz transformation R = (R%, i? 2 , R3Y, with Rj = —id Xj (—A)~^ 
being the jth Riesz transformation. 

Vp = RR-{u- Vu) - RR ■ (F. k ■ VF. fe ). 

In virtue of the boundedness of Riesz operator R in L p space for 1 < p < 00, we obtain that 

||Vp||a<C||Vtt||2||«||B + C||VF|| a ||J'||oo. 

For details about Riesz transformation see [TUJ [TT] . 
Thus from the first equation of (|1.3p we have 

n\\Au\\ 2 < ||«t||2 + ||« • V«|| 2 + ||Vp|| 2 + \\F. k • VF. fc || 2 

< |K|| a + |||Au|| 2 + C||u|| 2 ||Vu||* + C\\F\U\VF\\ 2 , 
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where the interpolation inequality ||«||oo < CIMI2 HAu^ nas been used. So we derive 

(3-5) ||Au|| a < C(||ut|| a + ||u|| a ||Vu||$ + ||F|U|VF|| a ). 

Next we derive the estimate of ||AF|| 2 . Applying A on the both sides of equation (|1.3|l and 
taking the L 2 inner product with Au and AF.fc respectively, we have 

(3.6) i^||A u ||2 + M ||AV M ||2 < \(A(u-Vu)-u-VAu,Au)\ + \(A(F. k -VF. k )-F. k -VAF. k ,Au)\, 

(3.7) iA||AF. fe ||| < |(A(« ■ VF k ) - u ■ VAF. k ,AF. k )\ + \(A{F k ■ Vu) - F k ■ VAu,AF k )\, 

where use has been of the facts 

(u ■ VAu, Au) = (u ■ VAF. k ,AF. k ) = 0, 
{F. k ■ VAF. kl Au) + {F. k • VA«, AF. k ) = 0. 
Next we estimate the right hand sides. By the communicator estimate (|2 . 10[) one has 
\{A{u ■ Vu) - u • VAu, Au)| < ||Au|| 2 [[A(u • Vu) - u ■ VAu\\ 2 < || Vw|[co|| Aw[[|, 
|(A(« • VF. fe ) - u ■ VAF. k ,AF. k )\ < ||AF|| 2 (||Vu|U|AF|| 2 + \\F\U\ VAu\\ 2 ) 

< UVulloollAf ||3 + CWFWUAFWl + |||VA«||1, 

\(A(F k • Vu) - F k • VAu, AF, k )\ < ||V«[[ <so ||AJ'||i + C||F|£,||AF||2 + |||VAu|£ 
For the second term on the right hand side of (|3.6p we estimate as follows 

\(A(F k ■ VF k ) - F k ■ VAF. k , Au)| < || Au|| 6 ||A(F. fe • VF. k ) - F. k ■ VAF fe || 6/5 , 

and 

\\A(F. k -VF. k )-F. k -VAF. k \\ 6/5 < ||VF|| 6 ||AF|| 3/2 < |[VF[[ 6 ||VF||f ||AF||f. 
So one has the estimate 

|(A(F fc • VF fc ) - F. k ■ VAF. fc , Au)| < |||VAu|| = + C\\VF\\* + C||VF||1||AF||§. 

Summing up (|3.6I) and (I3.7[) . and inserting the above estimates into the summation, we arrive at 
(3.8) 

_,_ u^. fc |||) + M ]| vAt*||i < v«Hoo + ll^lli, + 1] vj-iiiXH A«ni + II ap'HI) + C7]] v-F-nt 

We still have to estimate ||VF||6. Differentiating the second equation of (|1.3p with respect to Xi, 
one has 

dtdiF.k + d lU • VF. k + u • Vd l F. k = d l F. k ■ Vu + F. k ■ Vd lU . 
Multiplying both sides of the above equation by 6\diF. k \ 4 diF. k , and integrating both sides with 
respect to x over M 3 , it can be derived that 

(3-9) ^llVFll^CHVulUllVFllg + CllFlloollAullBllVFllg. 

Next we have to derive an estimate of ||Au||6- Using an argument similar to deriving the L 2 
estimate ||Ait||2 in (I3.5P we have 

Ml|Au|| 6 < \\d t u\\ 6 + llulUHVulle + C||^||oo|| Vi^|| 6 

( ' } ^llftVulla + CllulloollAulla+CHFlloollVFllB. 

Inserting estimates (13. 10)) to p.9[) one has 

(3 11) d7 l|VF|l « - C Q™°° + + \\d t Vu\\ 2 )\\VF\\t + CHFlUIMUIAuiyVFHi + C 

< CdlVulU + \\F\\l + \\d t Vu\\l + 1)||VF||* + CIIFH^IIullallAullJ + C. 
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Combining the estimates (13.81) and (I3.11[) we arrive at 
A(||A^ + ||AF. fc ||I + ||F|||) + M ||VA w ||l< 



CdlVulloo + IIVFII1 + \\FWl + \\d t Vu\\1 + l)(||Au||l + \\AF\\ 2 2 + \\F\\i) + C||F||i,N| 2 ||Au||? + C. 
Gronwall's inequality implies the H 2 estimates: 

||Au||! + ||AF. fc ||l + \\F\\i + nj WVAugds < exp |c(t) J (UVulU + \\d t Vu\\ 2 2 )dsj x 



(3.12) 



|Au(0)||l + ||Af. fc (0)||l + ||F(0)||* + C / (\\F\thh\\Au\\l + l)ds < co 



Based on the H 2 energy estimate the higher energy estimate can be obtained by bootstrap 
method as we did in section two. Thus the proof of the case n — 3 is completed. 
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